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Abstract 

In a recent paper it has been shown that for double extremal static spherical symmetric BPS 
black hole solutions in the STU model the well-known process of moduli stabilization at the horizon 
can be recast in a form of a distillation procedure of a three-qubit entangled state of GHZ-type. 
By studying the full flow in moduli space in this paper we investigate this distillation procedure in 
more detail. We introduce a three-qubit state with amplitudes depending on the conserved charges 
the warp factor, and the moduli. We show that for the recently discovered non-BPS solutions it 
is possible to see how the distillation procedure unfolds itself as we approach the horizon. For the 
non-BPS seed solutions at the asymptotically Minkowski region we are starting with a three-qubit 
state having seven nonequal nonvanishing amplitudes and finally at the horizon we get a GHZ 
state with merely four nonvanishing ones with equal magnitudes. The magnitude of the surviving 
nonvanishing amplitudes is proportional to the macroscopic black hole entropy. A systematic study 
of such attractor states shows that their properties reflect the structure of the fake superpotential. 
We also demonstrate that when starting with the very special values for the moduli corresponding 
to flat directions the uniform structure at the horizon deteriorates due to errors generalizing the 
usual bit flips acting on the qubits of the attractor states. 

PACS numbers: 11.25.Mj, 03.65.Ud, 03.67.Mn, 04.70.Dy 
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I. INTRODUCTION 



Recently striking multiple relations have been discovered between two seemingly unre- 
lated fields: the physics of black hole solutions in string theory and the theory of quantum 
entanglement within quantum information theory^"-. Further papers established a complete 
dictionary between a variety of phenomena on one side of the correspondence in the lan- 
guage of the other. This black hole qubit correspondence has repeatedly proved to be useful 
for obtaining additional insight into both of the two fields^"—. The main correspondence 
foundii^i^"-»i^ is between the macroscopic entropy formulas obtained for certain black hole 
solutions in supergravity theories and multiqubit and qutrit entanglement measures used in 
Quantum Information Theory. 

Apart from understanding black hole entropy in quantum information theoretic terms 
the desire for an entanglement based understanding for issues of dynamics also arose. In 
particular in the special case of the STU modeUi^i^ii^ it has been realized^ that for ex- 
tremal spherically symmetric BPS black hole solutions it is possible to rephrase the attractor 
mechanismi^ as a distillation procedure of entangled "states" of very special kind on the 
event horizon. Such states are of GHZ-type^^ or graph states^-ii^ well-known from quantum 
information theory. The basic tool for establishing this result was the introduction of a 
three-qubit state depending on the conserved charges and also on the moduli fields. 
Such "states" enjoy a number of remarkable properties^. The norm of this state having 8 
amlitudes is the black hole potentiaUi Vbh- The flat covariant derivatives with respect to 
the Kahler connection are acting on as bit flip errors on the qubits. At the horizon bit 
flip errors on are supressed for BPS solutions and for non-BPS ones they are not. The 
non-BPS solutions can be characterized by the number and types of bit-flip errors. 

However, in these investigations^'"^ establishing these results only double extremaP^ so- 
lutions have been considered for which the moduli fields are constant even away from the 
horizon. Since for this class of solutions is also constant clearly within the context of 
such solutions it is not possible to get any additional insight on the important question how 
the distillation procedure unfolds itself as we approach the horizon after taking the limit 
r — )■ with r being the radial coordinate. 

Luckily both in the BPS and non-BPS cases there exist more general static spherically 
symmetric solutions featuring the full radial flow in moduli space. For the BPS case these 
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are the well-known solutions based on harmonic functions-^ and for the non-BPS case sim- 
ilar results generalizing these ones have recently become available^i'^Si"— . In these works 
after solving the equations of motion one obtains the attractor flow z^{r) in moduli space. 
Hence employing the charge and moduli dependent multiqubit states^''' and using these so- 
lutions one might hope to get some additional insight into distillation issues by studying the 
corresponding flow |\E'(r)). 

The aim of the present paper is to investigate this distillation procedure in detail for 
the special case of extremal spherically symmetric black hole solutions in the STU model. 
We will use a special combination of the moduli fields, the warp factor and the conserved 
charges reminiscent of a 3-qubit state of Quantum Information Theory. We will call this 
creature a "three-qubit state" furnishing a representation space for the action of the duality 
group SL(2,M)^^ C Sp(8,M), though this terminology might be misleading. It will be 
obvious that our state has intimate connections to entities like the "fake superpotential"— 
and even possibly to the phase of the semiclassical wave function used in recent studies^^, 
however in this paper we will not elaborate on its physical meaning. Results on the origin 
of these 3-qubit states having some relevance on such interesting issues will be presented in 
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an acompanymg paper—. 

The organization of this paper is as follows. In Section II. we summarize the usual 
formalism of the STU model. In Section III. we introduce our moduli and charge- dependent 
3-qubit state, and recall results concerning the black hole qubit correspondence that we will 
need later. In Section IV. we reformulate the well-known findings concerning BPS solutions 
based on harmonic functions. Here we show that the "attractor at infinity"—'^ corresponds 
to a distillation procedure of a normalized GHZ state (dual to the usual one at the horizon^) 
at the asymptotically Minkowski region. In Section V. we study the flow |^(r)) for the 
non-BPS DO — D4 system answering the seed solution^i. Here we generalize further our 
3-qubit state by including also the warp factor into its definition. We show that the Fourier 
amplitudes of this state in the discrete Fourier (Hadamard) transformed basis satisfy a set 
of first order differential equations. Using the results of the previous sections in Section VI. 
for the non-BPS seed solutions we demonstrate how a GHZ state at the horizon emerges 
from a state characterizing the flow at the asymptotically Minkowski region. The attractor 
mechanism in this picture simply amounts to the fact that three amplitudes out of the 
seven nonequal nonvanishing ones are dying out as we approach the horizon. The remaining 
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amplitudes have the same magnitudes related to the macroscopic black hole entropy. The 
relative phase factors of these amplitudes are merely signs reflecting the structure of the 
fake superpotentiaUi>2i. In Section VII. using recent general results on the STU modelli we 
calculate the explicit form of the states at the horizon for both the available BPS and non- 
BPS solutions. Featuring the parametrizationiii^ revealing the fiat directions''^ we show 
that the role of the fiat directions in this picture is to deteriorate the uniform GHZ-like 
structure on the horizon. According to the results of our previous paper^ the differences 
between different types of solutions (BPS, non-BPS with vanishing^ and nonvanishing^ 
central charge) manifest themselves in applying bit flip errors to the relevant GHZ-like state 
of the BPS flow. In view of this result flat directions give a new twist on this picture, namely 
when starting the flow in one of the fiat directions the resulting state on the horizon will 
exhibit errors of more general type than the usual bit flip ones. Finally in Section VIII. we 
present our conclusions, with some calculational details left to an Appendix. 

II. STU BLACK HOLES 

In the following we consider ungauged N = 2 supergravity in = 4 coupled to n vector 
multiplets. The n = 3 case corresponds to the STU model. The bosonic part of the action 
(without hypermultiplets) (in units of Gat = 1) is 



Here J-"^, and *J^^, I = 0,l,2...n are two-forms associated to the field strengths J-"^^ of 
n + 1 U(l) gauge-fields and their duals. The z"" a = 1, . . .n are complex scalar (moduli) 
fields that can be regarded as local coordinates on a projective special Kahler manifold Ai. 
This manifold for the STU model is SL(2, ]R)/U(1)^^. In the following we will denote the 
three complex scalar fields as 




+ (ImA/}jJ^^ ■ J^^ + ReA/}jJ^^ ■ 



(1) 



2" = x'^ - ^^/^ 



a = 1,2,3. 



(2) 



With these definitions the metric and the connection on the scalar manifold are 



The metric above can be derived from the Kahler potential 
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as = dad^K. For the STU model the scalar dependent vector couplings ReA//j and 
ImA//j take the following form 
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We note that these vector couplings can be derived from the holomorphic prepotential 
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x^ = (x°,xV^ 



(7) 



via the standard procedure characterizing special Kahler geometry^. 

For the physical motivation of Eq. ([1]) we note that when type IIA string theory is 
compactified on a of the form x x one recovers N = 8 supergravity in = 4 
with 28 vectors and 70 scalars taking values in the symmetric space E7(7)/SU(8). This = 8 
model with an on shell U-duality symmetry E7(7) has a consistent N = 2 truncation with 4 
vectors and three complex scalars which is just the STU modeU^*^. The DO — D2 — DA — D6 
branes wrapping the various give rise to four electric and four magnetic charges defined 
as 

^ ' ^ 7 = 0,1,2,3 (8) 



^ Att Js2 ^ ' 



where 
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Qi=JfijF 

These charges can be organized into symplectic pairs 
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and have units of length. They are related to the dimensionless quantized charges by some 
dressing factors. Using the variable r = l/r = l/|x| and normalizing the asymptotic moduli 
as 

y^iO) = 1 x"(0) = B" (11) 

the dressing factors are essentially the masses of the underlying branes^. 

We are interested in static, spherically symmetric, extremal black hole solutions associ- 
ated to the (H]) action. The ansatz for the metric is 



ds2 = -e2^Mdt2 + e-2^(-)dx2 



(12) 



where the warp factor is a function of r = 1/r. Putting this ansatz into ([T]) we obtain a 
one dimensional effective Lagrangian for the radial evolution of the quantities U{t), -2"('r), 
as well as the electric and magnetic potentials^ 



and the constraint 
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+Gaa-^-; e^^VBHiz,z,P,Q) = 0. 

dr / dr dr 



(14) 



Here our quantity of central importance is the black hole potential Vbh which is depending 
on the moduli as well on the charges. Its explicit form is given by 
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where the matrices u = ReAf and fi = lmj\f are the ones of Eqs. ([5]) and ([6]). The explicit 
form of is 
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An alternative expression for Vbu can be given in terms of the central charge of = 2 
supergravity, i.e. the charge of the graviphoton. 
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(17) 



where for the STU model 

Z = e^l-'W = e^!\Q^ + z^Qi + z^Q^ + z^'Q^ + ^^z^^^po _ ^2^3pi _ ^i^3p2 _ ^i^2p3)^ (^g) 

and Da is the Kahler covariant derivative 

DaZ = {da + \daK)Z, (19) 

and W is the superpotential. 

Extremization of the effective Lagrangian Eq. f|T3|) with respect to the warp factor and 
the scalar fields yields the Euler-Lagrange equations 

ij = e^''VBH, z'' + Tlz'z'' = e^''d^VBH. (20) 

In these equations the dots denote derivatives with respect to r. These equations taken 
together with the constraint Eq. f[T^ determine the black hole solutions whose quantum 
information theoretic interpretation we are interested in. 



III. THREE-QUBIT STATES 

It is useful to reorganize the charges of the STU model into the 8 amplitudes of a three- 
qubit state 

|r)= ^ikj\lk3) \lk3) = \lh®\k)2®\3)i (21) 

Z,fcj=0,l 



where 




-P ) ^ ' ^ ' ^ \ _ \ ■'^000) Tool, Eoio, Eioo 
-Qq, Qi, Q2, Q3 I Vriii) riiO) Tioi, Foil 



(22) 



Notice that we have introduced the convention of labelling the qubits from the right to the 
left. Moreover, for convenience we have also included a factor ^ into our definition. The 
state |r) is a three-qubit state of a very special kind. First of all this state defined by the 
charges need not have to be normalized. Moreover, the amplitudes of this state are not 
complex numbers but real ones. As a next step we can define a new three-qubit state 
depending on the charges F and also on the moduli^'''. This new state will be a three-qubit 
state with 8 complex amplitudes. However, as we will see it is really a real three-qubit state, 
since it is U(2)^'^ equivalent to a one with 8 real amplitudes. 
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In order to motivate our definition of tlie new state we notice tliat^ 

Now we define the state as 

|*(- D) = e-r^ i^^ A « f^^^ A 8 -1 I ID. (24) 



Introducing the matrices 




a = 1,2, 3. (25) 



With this notation we have 

1^) = {S^®S2®Si)\V) = (U®U0l(){S3®S2®Si)\r). (26) 

This means that the state |\E') up to a phase for all values of the moduli is in the GL(2, C)^^ 

orbit of the charge state |r). Obviously the state |\E^) is an unnormalized three-qubit one 

with 8 complex amplitudes. However, it is not a genuine complex three-qubit state but 

rather a one which is U(2)^^ equivalent to a real one. This should not come as a surprise 

since the symmetry group associated with the STU model is not GL(2, C)^^ but rather 
SL(2,M)x3. 

Using these definitions we can write the black hole potential in the following nice form 

VBH = mf. (27) 

Here the norm is defined using the usual scalar product in ~ ® ® with complex 
conjugation in the first factor. Since the norm is invariant under U(2)^'^ our choice of the 
first unitary matrix of Eq. (125!) is not relevant in the structure of Vbh- We could have 
defined a new moduli dependent real state instead of the complex one |\E') by using merely 
the SL(2,M) matrices of Eq. (|25l) for their definition. However, we prefer the complex form 
of Eq. ( 126|) since it will be useful later. 

It is instructive to write out explicitly the amplitudes of our complex three-qubit state 

72^000 = e^/2^(^^^^^^), 72^111 = -e^/2y^(z^z^z^), (28) 
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with the remaining amphtudes arising by cychc permutation. Notice also that we have the 
property (reahty condition) 

^000 = -"^iii, ^110 = -"^ooi, ^101 = -"^010, ^011 = -^100, (30) 
which can be written as 

1^) + (ai ® ai ® ai)|^) = (31) 
via the bit flip operator (Xi. Using this in Eq. (!27|) we can write Vbh in the alternative forn>2^ 



Vrh = e 



K 



\W{z\z\z')\^ +\W{z\z\z')\^ + \W{z\z\z^)f + \W{Y',z^z^)f^ . (32) 



As a motivation for our particular definition for we notice that 

^2*000 = Z, (33) 

i.e. the amplitudes and \l/ooo are related to the central charge and its complex conjugate. 
The remaining amplitudes are simply arising by conjugating one or two moduli. Moreover, 
one can show^ that the flat covariant derivatives with respect to the moduli act on our 
as bit flip errors. Explicitly we have 

D^l^) = (/®(T+®/)|^), D^l^) = (/®a_®/)|^), (34) 
D^\^) = (a+®/®/)|^), Z^^l^) = ((T_ ®/®/)|^). 
Here the operators a± act as 

a+|0) = |l), fT+|l) = 0, a_|0)=0, fT_|l) = |0), (35) 

and the flat covariant derivatives are defined as = —2iy^Di, Dj = 2iy^Dj e.t.c. where 



D,W{z\z\z') = ^^_!''f'f\ DjW{z\z',z') = e.t.c. (36) 



Hence the flat covariant derivatives are acting on our three-qubit state as the operators 

of projective errors known from the theory of quantum error correction. Alternatively one 
can look at the action of the combination Dq, + -D^ 

(Di + D^)|^) = (/®/®(Ti)|^), e.t.c. (37) 



According to our previous paper- it is illuminating to use the discrete Fourier transform 
of our three-qubit state At the event horizon the moduli are stabilized due to the 

attractor mechanism, and their stabilized values can be expressed in terms of the charges. 
These stabilized values give rise to entangled states on the horizon of very special form. 
It was shown'^'^ that for BPS solutions these states or of generalized Greenberger-Horne- 
Zeilinger (GHZ) fornt^^, and for the simple non-BPS solutions^^ they are graph-states^ 
well-known from quantum information theory. 

The discrete Fourier (Hadamard) transformation is implemented by acting on by 
H ® H ® H where 

H=^\ ^ \ (38) 



72 



1 -1 



Hence the Fourier transformed basis states are defined as 

|0) = ^(|0) + |1)) = i/|0), |i) ^ i=(|0) - |1)) = H\l). (39) 

Since HaiH = and Ha^H = ai the operator ai is acting on the Hadamard transformed 
base as a phase (sign) flip operator and vice versa. The important corollary of this obser- 
vation is that in the theory of quantum error correction once we have found a means for 
correcting bit flip errors using a discrete Fourier transform the same technique can be used 
for correcting phase flip ones. It was shown^ that the phase flip errors in the discrete Fourier 
transformed base correspond to flipping the sign of certain charges of the symplectic vector 
F. For BPS solutions these phase flips are supressed, and for the simple non-BPS solutions 
only errors of very special kind are allowed. 

Now we would like to gain some more insight into these interesting results by studying 
the solutions even away from the horizon. For this the full solution of the flow in moduli 
space is needed. For BPS solutions we can use the well-known results^*^ and for the non- 
BPS solutions the recently found seed solutions^S^^ and the most general non-BPS flows^ 
generalizing the simple non-BPS solutions^Si^^i^. 

Our main calculational tool will be to consider the discrete Fourier transform of our |\I') 

1^) = {H®H®H)\^!) = {V®V®V){S3®S2®Si)\r) (40) 

where 

V=\ \ (41) 
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is just i times the usual phase gate known from quantum information theory. Hence from 
Eq. fHOl) we see that the Fourier transformed state is up to some important complex phase 
factors is lying on the SL(2, R)^^ orbit of the charge state |r). Now solving Eqs. (^U^ with the 
constraint of Eq. f[T^ yields the flow z"-{t) on moduli space. To this flow we can associate 
a corresponding one |\1'(t)) of the Fourier transformed state. Our aim for the following 
sections is to look at the structure of |^E'(t)) for BPS and the non-BPS seed solutions. 

IV. BPS SOLUTIONS 

In this section we would like to study the behavior of BPS solutions in the three-qubit 
picture. In particular we would like to see how the three-qubit state of Eq. (12^ behaves 
as a function of r, answering to the flow z"'{r) in moduli space. This three-qubit picture 
is natural as the [/-duality group is SL(2,]R)^^ which is a subgroup of Sp(8,M), hence we 
expect that the usual symplectic invariants occurring in the formalism of the STU model 
should boil down to the corresponding SL(2,R)^^ i.e. three-qubit ones. 

As a first step in order to present the solutions capable of incorporating a wider range of 
asymptotic data with i?-fields we define a set of harmonic functions as 

'H{T)=T + Tr, i.e. = p' + P'r, Hj = Qj + Qjt /, J = 0,1, 2, 3. (42) 

We can alternatively encode the asymptotic data ('H(O) = F) into a three-qubit state 

irr (^000, Tool) Toio, Fioo\ 1 ( P ^ P , P , P \ 

|r)= 2^ ^ikj\lkj), I _ _ _ \ = _ _ _ _ \ . (43) 

i,k,j=o,i yriii, Fiio, Fioi, Foiiy v2 y— Qq, Qi, Q21 Qzj 

which plays a role similar to the charge state of Eq. f l2T]) . Notice that unlike for F in Eq. ( H2]) 
in the state |F) we also included a factor of ^ for convenience. 
For later use we also define a r dependent three-qubit state as 

\n{r))^\V)^r\V). (44) 

This state in the limits r — > and r — > 00 characterizes the asymptotic data and the charge 
configuration respectively. 

Let us now define Cayley's hyperdeterminant^>2^ D(\il})') for an arbitrary three-qubit 
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state IV') = ^ikj=oii^ikj\^kj) with amplitudes 

^000, V^OOl, ^010, V^lOO j ^ j ^0, ^1, ^2, ^"4 . ^^^^ 
^111, V'llO, ^101, ^011/ \i^7, A, i^b, 



as 



-2[(V'1^6)(V'2^5) + i'ipMii^M + (V'2^5)(^4V'3)] + 4^oV'lV'2V'4 + 4V'7V'6^5^3- (46) 

D{\iIj)) is permutation and SL(2, C)^^ invariant and under the group GL(2, C)^^ transforms 
as 

D{^) ^ (DetG3)2(DetG2)^(DetGi)2D(V^), Gg ® G2 ® Gi G GL(2, C)^^ (47) 

Notice that for a general charge state |r) such as the one of Eq. fl2Tl) we have 

- ADi\T)) = hiV) (48) 

where /4(r) is the usual quartic invariant known from studies concerning the STU model 

h(T) = - {P'Qir + i[{P'Q^){P'Q2) + {P'Q2){P^Q,) + {P'Qi){P'Q,)] 

+ AQoP'P^P^ - ^P^QiQ2Q3. (49) 

According to the general theory the data giving rise to BPS black hole solutions are 
incorporated into a 'H(r) subject to two constraintai^i^ 

/4(r) = -4D(|r)) = 1, (r, r) = p'q, - Qjp' = o. (so) 

As far as these constraints hold we can completely characterize any /4(r) > solution by 
generalizing the attractor equations to the so called stabilization equationsi^. The warp 
factor of Eq. ([12]) is 

^-.uir) ^ ^ -m\H{r)))- (51) 

Since for normalized states the quantity 

0<ri23=4|Z}(|^))|<l (52) 

is characterizing the tripartite entanglement of three-qubit systems^^ we observe that the 
warp factor can be regarded as a r = dependent entanglement measure describing the 
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tripartite entanglement of a state encapsulating the details of the charge configuration and 
the asymptotic values for the moduli, i.e. 

e-^^(^) = ri,3(|H(r))). (53) 

In this picture the first constraint of Eq. fl50|) means that for the state of Eq. (jUj) the 
asymptotic value of the "three-tangle" is normalized to one. 

In order to find an entanglement based meaning for the second constraint of Eq. fj50]) 
notice that the 16 quantities of the states |r) and |r) can be organized into a four-qubit 
state. Indeed let us define a state I7) with its 16 amplitudes •jmikj given by 

lOlkj = ^Ikj, lllkj = ^Ikj- (54) 

then a four-qubit entanglement measure invariant under SL(2,C)^'^ is given by-^ |o"i234| 
where 

cri234(|7)) = 7o7i5 - 7i7i4 - 72713 + 73712 - 74711 + 75710 + 7679 - 7778 (55) 

where for simplicity we again used the decimal labelling. Now using the definitions as given 
by Eqs. fl2T]) and (H3|) it is straightforward to check that 

(r,r) = ai234(|7)) (56) 

meaning that our second constraint is equivalent to the vanishing of the entanglement in- 
variant (T1234. Hence we conclude that both of the constraints describing the BPS black hole 
solutions have a characteristic meaning in our entanglement based reformulation. Recalling 
also that the value of Ti23(|r)) is related to the the entropy of the BPS STU black hole^i^*^ 
we can summarize these results as 



^ = vrv/ri23(|r)), ri23(|r)) = l, ai234(|7)) = 0. (57) 

Now in order to set the stage for the following generalizations we review and slightly 
extend the known results concerning the distillation procedure for BPS solutions^. In the 
following for simplicity we consider the DO — DA system. In this case we have Qq > and 
P* > but = Qi = 0. Generalizing the simple BPS solution we also include non-trivial 
B-fields as foUows^i. First we define our state |r) with amplitudes for Eq. (H3|) as 

P° = ^ sin 5, P''''' = ^ (cos 6 + sin 6B^'^'^), (58) 

13 



Qi = ^(sin(5[l - B'^B'^] - cos5[i?^ + -B^]), and cyclic permutations (59) 
v2 

= ^ [(5^ + 52 + 5^ - B^B'^B^) sin 5 + (1 - 5^5^ - 5^5^ - 5^5^) cos 5] . (60) 
Here 5^, B^, B^ are related to the asymptotic values for the moduli as 

zX=,^B''-i, (61) 

i.e. the asymptotic volume moduli are normalized, but we keep the asymptotic i?-fields as 
free variables. 

In terms of our three-qubit states there is a nice way of understanding this choice for 
and Qj and also the meaning of the additional parameter 5. First just like in Eq. fl26|) we 
define a new moduli dependent state \'^{t)) as 

|^(r)) = (53(r) 52(r) 5i(r))|r), (62) 



where for the definition of 5i,2,3 see Eq. (1251) . We would like to see how this state behaves 
at the asymptotically Minkowski region. A straightforward calculation using Eq. ( 16T|) shows 
that 

1^(0)) = ^(e-*^|000) - e^^llll)), (63) 
v2 

i.e. the parameter 5 is related to the phase of a generalized GHZ state. Notice that thanks 
to our inclusion of the factor ^ in the (1431) definition of |r) this state is normalized. On 
the other hand, according to Eq. ( 1271) a similar inclusion for the definition of the companion 
state |r) (Eq. ( l22l) ) has also fixed the norm of the state |\E') to be the black hole potential. 
Knowing that for the normalized generalized GHZ state of Eq. (163 p ri23 = 4|D| = 1 by virtue 
of Eqs. P7|) . (H9|) and (!62|) we immediately get /4(r) = 1, i.e. the first of our constraint is 
satisfied. 

In our recent paper- it was shown that the attractor mechanism for STU BPS black holes 
can be reinterpreted as a distillation mechanism of a GHZ state at the horizon (r = 00) of 
the form 

\^>{oo)) = (/,(r))i/^i=(e-^^|000) -e^^llll)), (64) 
where the explicit expression of A is given by- 
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where for the definition of /4(r) see Eq. (H9|) . Comparing Eqs. (l63ll and (IMIl we see that we 
have an attractor at the horizon (r = oo) and an attractor at the asymptotically Minkowski 
region (r = 0). These attractors can be described by the fiows of the charge and moduli 
dependent states |^(t)) and |^(t)) respectively, each of them producing a maximally en- 
tangled GHZ state at the attractor points defined by r = oo and r = 0. The attractor at 
r = giving rise to the GHZ-state of Eq. (163|) is the well-known attractor at infinity^i^ a 
map between the 6 real moduli and the 8 constants in the harmonic functions subject to 2 
constraints. We can summarize these considerations by noting that for the relevant fiows 
and attractor points we have 

||^(0)|p = 1, ||^(oo)||2 = Vbh{oo) = ^/hif). (66) 

Now we compare the fiows |^E'(t)) and |^E'(t)) by calculating 

(^(r)|a3®a3®a3|^(r)). (67) 

This quantity is the transition amplitude between the phase- fiipped state {a^®a^®a^)\'^{T)) 
and the one \ '^{t)). Using Eqs. (J26ll and (1621) and the fact that Wa^U = —02 and S'^ aiS = a2 
for S e SL(2,M), we get 

(^(r)|(T3®(T3®(T3|^(r)) = -(r|(T2 ®CT2 ® (T2|r) =2(1,1) = 2(Ti234(|7)) =0. (68) 

This shows that the phase (sign)-fiipped version of ^I'(t) is always orthogonal to the com- 
panion state ^(r). It is also clear that this amplitude is of purely topological in origin. 
(E.g. in the type JIB duality frame it is related to the intersection product^^ on T^.) Due 
to the fact that Ha^H = ai, where ai is the bit fiip operator we can alternatively conclude 
that the bit fiipped version of the Fourier transform of one of the states is orthogonal to the 
Fourier transform of the other. Recall also that according to Eq. (1371) such bit fiip errors are 
related to the action of the fiat covariant derivatives with respect to the moduli. 

As a last application of the vanishing of the amplitude of Eq. ( I68l) . Eqs. (1631) and ( l33ll 
give the meaning of the parameter 6 as the phase of the central charge Z (up to a shift 
by vr). These results shed some light on the quantum information theoretic meaning of the 
second constraint of Eq. (!50|) . 

In closing this section we present the well-known solution of the stabilization equations 
giving the moduli fields as a function of r— "^"^ 

1 ^ -HiH^ + H,H^ + H2H2 + H:,H^ - le-''^ 
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with the remainig equations for and are arising after cychc permutation of the numbers 



given by Eq. ( ISTl) . It is straightforward to check that the solution z°-{t) satisfies the correct 
asymptotic behavior of Eq. ( l6Ti) and 



By virtue of Eqs. (!28|) - (!30|) and (!6T]) it is clear that the asymptotic values of the 8 amplitudes 
of 1^(0)) are generally non zero. However, the flow z°'{t) in moduli space giving rise to the 
flow |^(t)) results in a GHZ state at the horizon having merely 2 nonvanishing amplitudes 
(see Eq. (El])). 

This distillation process is the one that was studied for BPS solutions^. Later a non- 
BPS generalization was also given^. In this case the distillation process gives rise to graph 
states^^ at the horizon. However, the analysis in these papers was restricted merely to double 
extremal solutions^^ for which the moduli are constant even away from the horizon. Hence 
in these studies the important question of how the distillation process becomes unfolded as 
r changes have not been addressed. Our aim in the next section is to investigate such issues 
by studying the non-BPS flow explicitly. 

V. THE FLOW FOR THE NON-BPS D0-D4 SYSTEM 

In this section we study the flow |^(t)) answering the full radial flow z°'{t) obtained 
for the 5 parameter family of non-BPS seed solutions of a — DA system discribed by 
Gimon et.al.— More precisely it turns out to be rewarding to study the properties of a 
related flow |x(t)) instead by multiplying |^'(t)) by the warp factor. Then we show that 
the Fourier amplitudes of this new state satisfy a first order system of differential equations 
in accordance with our expectation coming from previous studies^ii^^iSl. 

First we address the quantum information theoretic aspects of the seed solution. For the 
DO — DA system we chose Qo < and P'* > 0, a = 1, 2, 3. Let us define 



1, 2, 3. Here = + PV, Hj = Qj + Qjt, with /, J = 0, 1, 2, 3, 



and the warp factor is 




(70) 




(71) 



From Eq. (120|) the equations to be solved are 



A (4e-2*») = 2e' 
dr ^ ' 



OVbh 

dta 



a = 1,2,3 



(72) 
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4>a+{iae-''^ Y = 2e'^^, a = 1,2, 3, (73) 

U = e'^VsH- (74) 
Moreover, according to Eq. f|T4l) we also have the constraint 

U' + \Y.{^" + [iae-^^r)=e'^VBH. (75) 

a 

As a first step we define a new three-qubit state by incorporating also the warp factor as 

|x(r)) = e^W|vl/(r)) = e^^^^^s ® 5^ ® 5i)|r). (76) 

This state depends on the charges, the moduli and the warp factor. The discrete Fourier 
transform of this state (for the definitions see Eqs. fl25l) . fl38|) . fHOj) and f HT]) ) 



|x(r)) = {H®H^ H)\x{t)) = e^^^\V ®V ® V){S^ ®S2® 5i)|r), (77) 

will play a particularly important role in the following. 

For later use we write out explictly the amplitudes of this Fourier transformed state | %(''")) 

Xooi = \\h\'^e^^^^^^\ xoio = \\h\^eP^^^^^\ Xioo = ^|/4|^e^+^^+'^^ (78) 

%io = ^|/4|^e^+<^ni2 + ^3), %oi = ^|/4|^e^+^Hti + t3), ^Xoii = i|/4|^e^+<^^(ti +t2), 

(79) 

Xooo = 0, Xiii = -2l^4|^e'^(l + tit2 + t2t3 + t3ti). (80) 

Here 

/3 = - i(0i + 02 + 02), I/4I = -h = -AQoP'P^P' > 0. (81) 
Now in terms of these amplitudes the equations to be solved can be written as 

- {i^^"^^) = -mYalx), (82) 

4>a+{iae-'^^Y = 2{x\Za\x), (83) 
U={X\X), (84) 

where the operators Zi = I ® I ® a^, Z2 = I ® ® I and Z^ = ® I ® I describe phase 
fiips of the ath qubit in the Fourier transformed base. The quantities Ya are defined using 
a2 accordingly. Notice also that the operators \{1 + Za) where 1 is the 8x8 unit matrix, are 
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projection operators. Hence by adding the half of Eqs. (l83i) to Eq. (18^ gives rise to three 
equations which contain merely three amplitudes (xo = 0) on the right hand side. These 
manipulations also justify the introduction of a new variable^Si^i 



aa = U + ^(j)a, 

Using this new variable instead of Eqs. ( 183|) we can use the ones 



(85) 



"2-21 -t2e 



"3 - 2 ( -he 



-4>2 



-4>3 



^ixl + xl-xl): 



nxl + xl~xl) 



(86) 

(87) 
(88) 



where from now on we use decimal labelling for our amplitudes. For Eq. fl82|) we have the 
form 



/a. 

dr 



v 







^ 

dr 








_d_ 

dr 



V 



\2 



X7 -X4 -X2 

Xi Xi -Xi 
X2 -Xi X-j 



X5 

\Xzj 



(89) 



and the sum of Eq. (18^ and —\ the sum of Eqs. (183|) gives 



/3 + 2$^(^tae 



4X7-2(X^ + X5 + X3)- 



(90) 



Now, the crucial observation which enables an explicit construction of the seed solutions is 
the fact that in the Fourier transformed basis we have seven nonvanishing amplitudes, and 
the constraint Eq. (1751) is consisting of squares of seven terms. Using the decimal labelling 
for the amplitudes of |x) this constraint can be written in the form 



/j' + ^ {aa + dfe - dc)' + J](e-'^na)' = 4(x? + X? + X2 + X4 - Xe - Xs - x\)- 
A natural choice to satisfy this constraint up to a sign is 

± X7 = ±X6 = ^e-'^^ii, ±X5 = ^e-<^^i2, ±X3 = ^e-^^tg, 

± Xi = ^(«i - «2 - «3), ±X2 = ^("2 - ai - as), ±X4 = ^(as - "1 - a2)- 



(91) 



(92) 

(93) 
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Comparing these with the exphcit form of the Fourier amphtudes of Eqs. ( !78|) - (l80|) with 
Xi,2,4 rewritten as 



_ _| r|l/4 02+03-01 
Al — 2 I 4| ; 



X2 



2 2 



we get the following set of first order differential equations 



4 = T|/4r/'(4 + tc)e'""-"^-"^ 



(95) 

(96) 
(97) 



where a ^ h ^ c. The solutions to these equations with the upper sign choice were given 
in Gimon et.al.— Before recalling these solutions we show that these choices automatically 
solve Eqs. (I86|) - (p0|) . Let us substitute instead of ^iaC''^" occurring in these equations the 
Fourier amplitudes of Eqs. f l92|) . Then we get 



ai = 2{xl + xl), a2 = 2{xl + xl), a3 = 2{xl + xl)- 
In Eq. fl89|) ve also replace x? by |/3 (we have chosen the upper sign) to get 

/i, 

2 

\ -X2 



(98) 
(99) 



d 

d7 



X5 
\X3j 



2(/3 + 0l) -X4 -X2 

-Xa |(/3 + 02) -Xi 

-Xi 10 + <P3)J 



X5 
\X3j 



(100) 



Now using the explicit expressions for the Fourier amplitudes it is easy to check that these 
equations are indeed satisfied. 

An interesting possibility is to write down these equations as first order equations for the 
Fourier amplitudes. From Eqs. f l98l) and (11001) we get 



(101) 





fxr\ 




A/3 


-X6 


-X5 


-X3 


\ 




d 


X6 


= 2 





|(/3 + 0i) 


-Xi 


-X2 




X6 


d7 


X5 







-X4 


i(/3 + 02) 


-Xi 




X5 










-X2 


-Xi 




>3)/ 
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Similarly from Eq. ( p9|) using Eq. (IMj) the corresponding equation is 

/ 

2 



d 

d7 



X2 
VX4/ 



v 



-Xi 
-X2 
-X4/ 



X2 
VX4/ 



(102) 



An alternative form for Eq. f llOip can be given by noticing that 

/3 + 0, = 2f/ - (dfe + dfe - da) =2U+ 1/41^/^6""+"=-"" = 2f7 + 2X1,2,4, 

hence we have 



(103) 









d 








= 2 


d7 


X5 





/v 



X7 



-X6 -X5 



-X3 \ 



X6 
\X3/ 



(104) 



f/ + Xi -X4 -X2 

-X4 + X2 -Xi 

.0 -X2 -Xi u + X4/ 

Eqs. (11021) and (I104p are first order equations containing the Fourier amplitudes of our 
charge, moduli and warp factor dependent states and U. However, we can express U in 
terms of some of the x s as follows. Let us define 



^ = \{X7 - Xi ~ X2 - Xi) 



(105) 



Then from Eqs. (fT02|) and (fTOD we get 

= ixi)' + iX2)' + (X4)' + ixr)' - - (X5)' - ixs) 

a T 

Hence fom Eq. ( 184]) 

The new quantity W 
W 



— U = w = e^W. 
dr 



^1 - ^2 - ^4) 



^(^111 - ^ 



001 



010 



100 J 



(106) 



(107) 



(108) 



is the fake superpotentiaUii^i where it is easy to check that its explicit form coincides with 
the negative of the one as given by Eq. (6.8) of the paper by Bellucci et.al.^^ Notice also 
that the fake superpotential contains only Fourier amplitudes of odd parity of the charge 
and moduli dependent three-qubit state This will be of some importance in the next 
section. 



20 



Since we have U 



\{X7 - Xi 



X2 — Xa) = w we can write Eq. fll04p in the final form 



d7 







X6 






= 2 


Xb 





^ X7 -Xd -X5 -X3 \ / X7\ 



x& 



(109) 



w + xi -Xi -X2 
-X4 W + X2 -Xi 
VO -X2 -Xi w + XiJ 

Hence Eqs. (11021) and (11091) show that in the case of the seed solution for the non-BPS Z ^ Q 
DO — DA system the r derivatives of the Fourier amplitudes of our charge, moduli and warp 
/actor-dependent state |x) can be expressed entirely in terms of the Fourier amplitudes. 



VI. THE ATTRACTOR MECHANISM AS A DISTILLATION PROCEDURE 

In this section we would like to demonstrate how the radial fiow studied in the previous 
section gives rise to the distillation of a special three-qubit state at the black hole horizon. In 
order to see this procedure unfolding all we have to do is to use the solutions of the first order 
equations Eqs. fl95l) - fl97j) to obtain explicit expressions for the Fourier amplitudes Xifcj (''")• 
It means that starting from the asymptotic values X«fcj(0) in the region with Minkowski 
geometry at the limit r — )• oo we obtain the ones xikj{oo) at the horizon with adS2 x S"^ 
geometry. The solutions of Eqs. ( |95l) - (l97j) are^ 

6°-+""-"- = , ,1^, ,,,, = ^, ay^b^c, a, 6, c= 1,2, 3, (110) 

ta = ^, = -ho-.rS^, ho = -do - Ihl'^^T (111) 

hbK hih2hs 

where 

17" |V4 \T |l/4 

where for the definition of I/4I see Eq. ( ISTl) . Notice that according to Eq. (fTTl) and ( |7T]) for 

this 5 parameter solution we have 

R _ ,>-2C/ .... 
^"^ = ^--1 TT' e-^"" = -hoh,h2h-B^ (113) 

hence B = x^{0) = x^{0) = x^{0). 

In order to write down the explicit form of the amplitudes of our charge, moduli and 
warp factor dependent Fourier transformed state it is useful to introduce the new harmonic 
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functions 
and the warp factor 



^ + ^'V = p^^„(r), H,{T)^-^il + B') + QoT^-j^^ho, (114) 



g-4[/(r) ^ -AHo{t)H\t)H\t)H\t) - B\ (115) 



Using these results we get the following results for the Fourier amplitudes for our charge, 
moduli and warp factor dependent state \x) 

pi p2 d3 



Xooo(t) = 0, Xooi('r) 



2m {tY 



Xoio('r) 



Xiio(r) = 


_lp2Lr(r) 

2 


Xioi(r) = 


_le2l/(r) 

2 


Xoii(r) = 


2 



2ii-2(^) ' 

p2 p3 



Xioo(t) = 



2/^3 (r) 



pi p3 



pi p2 



Xin(r) = ^e^^(^) 



4goi^^(r)i/^(r)i/3(r)-B^^ — 



po 



a=l 



mir) 



, (116) 

(117) 
(118) 
(119) 

(120) 



From this the components of the charge and moduli dependent Fourier transformed state 
|^(t)) are 

^,,,(r) = e-^«X^.,(r). 



Since 

1 

lkj=0 

one can show that in the asymptotically Minkowski region we have 

lim|^(r)) = -^(^P^IOOI) + P'^ \010) + \100) - iB{P^ + P^)\110) 

- iB{P^ + P^)\101) - iB{P^ + P^)\011) 

+ [Qo-B'(p^ + p^ + p^)]\in)y 

On the other hand at the horizon we have 

lim |^(r)) = |/4|^/4f|001) + I010) + |100) - Im))- 



(121) 



(122) 



(123) 



(124) 
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This result shows that if we "start" asymptotically with the state of Eq. fll23p with seven 
nonvanishing generically different amplitudes we end up with the state of Eq. ( I124p having 
merely four nonvanishing amplitudes that are the same up to a sign. Notice also that the 
four nonvanishing amplitudes are having states with odd parity. This reminds us of the 
structure of the fake superpotential Eq. fllOSp . This is as it should be since we know^^ that 
in the near horizon limit should give the square root of —I 4, > 0. Indeed, since we have 

lim ^i,2,4(r) = - lim ^r^r) = hh\'/\ (125) 

lim^^oo VV^(t) = 1/4!^''^ = \/ —AQqP'^P'^P'^ which is the correct value. According to 
Eqs. ( I48l) and ( l52l) in our three-qubit interpretation 

W2(oo) = y^ri23(|^(oo))), (126) 

i.e. the square of the fake superpotential on the horizon is just the entanglement measure of 
the state of Eq. fll24p . It is important to realize that the components of the fake superpo- 
tential are precisely those amplitudes of our 3-qubit state which are not dying out as we are 
approaching the horizon. According to our previous results this also works for the BPS case 
(and as will be seen in the next section even for the non-BPS case with vanishing central 
charge). In order to verify this just recall that for BPS solutions \Z\ plays the role of the 
superpotential which is according to Eqs. fl28|) and ( j64|) is again related to the amplitude 
which is not dying out in the attractor limit. 

VII. ENTANGLED STATES OF GHZ TYPE ON THE HORIZON 

In the previous sections we have studied the distillation process in the special case of the 
non-BPS Z 7^ seed solution. Clearly similar results can be obtained for the most general 
non-BPS solutions^^ and the non-BPS ones with vanishing central charge^. In this section 
however, our main concern will be to present the explicit forms of our 3-qubit state |^'(t)) 
at the horizon. Of course the states we expect to show up are again GHZ-like states, but 
the new subtlety worth investigating in this context is the appearance of flat directions^. 
As we have mentioned we have to make distinction between three different cases |-BPS, 
non-BPS Z = 0, and non-BPS Z 7^ solutions. 

In this section we use the , qi quantized charges instead of the P^, Qi dressed ones. 
This is because we would like to use the most general non-BPS Z 7^ solutionli which 

23 



has been produced by using an U-duality transformation acting nicely on such quantized 
charges. The dressed charges are rescaled quantities related to the quantized (undressed) 
ones via factors coming from the asymptotic volume moduli. Calculating the moduli z"" using 
the f[T5]) black hole potential with p^, g/, the asymptotic volume of the tori are nontrivial 
y''(0) = Va- In order to get ?/'^(0) = 1 as in Eq. (II ip we have to rescale the charges with 
real positive dressing factors. (For the definitions of these factors see the paper of Gimon 
et.al.—) The dictionary between the two conventions, i.e. , qi with 5" = z"'Va and , Qi 
with is then effected by the correspondence 

Vbh{P'. Qi, z") = GnVbh{p', Qi, z") (127) 

where is the D = A Newton constant. The quartic invariant defined in Eq. fj49|) can also 
be written in terms of the quantized charges, hence 

h{P',Qi)=G%h{p',qi). (128) 

In this section we denote I4 = I^ip^ , qi)- One can check that not only the norm of \E' i.e. the 
square root of the black hole potential, but also our three-qubit state scales simply with 

\^{P',Qj,z^)) = y^\^{p',qj,r))- (129) 



A. The BPS case 



After these techniqualities first we turn once again to the BPS solutions. The black hole 

charge configurations supporting the |-BPS attractors at the event horizon are the ones 
satisfying the following set of constraintsii 

h > 0, pV - p^qc > 0. (130) 

In this case the general |-BPS attractor flow solution is^iii^ 

exp{-4U) = h{h^,hi), (131) 

x-(r) - - (132) 
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Here the indices a, b, c are distinct elements of the set {1, 2, 3}, and no summation is imphed 
on them, on the other hand for the indices I = 0,1,2,3 summation is understood. The 
undressed harmonic functions are defined similarly to the f l42p dressed ones 

h'{r) =f+ p't, hrir) =qj + qiT. (134) 

In the horizon-limit we have 



lim x" 

T— >00 



lim y" 



(135) 



2(pfepc _ pOg^^ ' » 2(pV - P°9a) ' 

In order to obtain |\^(t)) on the horizon, we have to apply (P (g) P (g) P){Ss Cg) 5*2 (X) 5'i) on 
the charge vector and taking the limit r — )■ oo. For notational simplicity we introduce the 
abbreviation 

1^) = lim |*(r)). (136) 

T— >00 

Then using the identity 

4(pV - p°qi){p^p^ - p^q2){p^P^ - /gs) = {p^fh + (2pVV - P^P^Qi f (137) 

we get 



1^) 



—I 



-'4 



a(|000) + |110) + |101) + |011))+2/3(|lll) + |001) + |010) + |100)) . (138) 



2 + 

The discrete Fourier transform of this state is 



(3 — ia 



:|000) 



[3 + ia 



a" 



Hill) 



(139) 



with 



a = ^/\h\p°, 13 = 2pVV - pVg/, (140) 

in accordance with our previous results^ and Eqs. and (1651) . It is important to realize 
at this point that our state at the horizon can alternatively be written as 



-iarg(Z) 



|000) -e*^''§(^)|lll) 



(141) 



where the quantities \Z\ and arg(Z) now refer to the magnitude and phase of the central 
charge at the horizon. Recall that arg(Z(r)) satisfies the following equation 



^arg(Z(r))+^(r) = 0, 



2^7^' 

a=l ^ 



(142) 



i.e. A is the Kahler connection. Hence the relative phase factors that show up in the GHZ 
states of Eq. and (I139P are just the attr actor values for the phase of the central charge 
governed by Eq. ( I142p . 
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B. The non-BPS Z = case 



The non-BPS Z = solutions^iii can be obtained from the |-BPS ones simply by changing 
the sign of any two imaginary parts of the moduh. This yields the change of the f ll30p |-BPS 
constraints 



h > 0, p^" - p^Qc > 0, pY - P Qa < 0, pY - P Qb < 0. (143) 

During the calculation of |\&) the moduli appear only in the Sa matrices. We can carry out 
the sign flip of some y"' with the Pauli matrix: 



1 



-y'' 



X-x" 1 



(144) 



Because of this |^'23) = (erg (g) as (g) /)|^) and 1^23) = {H H H)\-^23) = (cxi (g) di (g) J)|*), 
where the indices of |\E'23) denote which moduli have been conjugated. By virtue of these 
observations we get 



1^12) 



1^23) 



1^13) 



-ill 



a(|000) - |110) - |101) + |011)) + z/3(|lll) - |001) - |010) + |100)) 

(145) 

a(|000) + |110) - |101) - |011)) + i/3(|lll) + |001) - |010) - |100)) 

(146) 

a(|000) - |110) + |101) - |011)) + i/3(|lll) - |001) + |010) - |100)) 

(147) 



The discrete Fourier transform of these states is 

|011) 



1^23) = n 



1 1 



V2 



/3- 


- ia 




+ a2 




- ia 




+ a2 


/3- 


- ia 




+ a2 



1110) 



13 + ia 
(3 + ia 

7WT~ 



:|100) 



:|001) 



a^ 



B + ia , 
|101) |010) 



(148) 
(149) 
(150) 



Note that a and (3 are the same as in the |-BPS case however, now the charge configuration 
should be compatible with the restrictions of Eq. (I143p . 
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We can also write these states in the form of Eq. (11410 . For example singleing out the 
first qubit we get 

1^23) = |^i|^[e^'^''s(^i)|110) -e-^^''s(^^)|001)], (151) 

where 

Za = DaZ = -2iy''DaZ, (152) 

for the definition of DaZ see Eq. (fT9|) . Here as in Eq. (1141 p the quantities \Zi\ and arg(Zi) 
refer to their attractor values. Generally these quantities are r dependent. For example 
arg(Zi(r)) satisfies the following equation 

i^arg(Z,W) + AM = 0. M^U% -%-%). (153) 
dr 2 V r / 

Hence just like the phase of the central charge arg(Z) for |-BPS solutions for this non-BPS 
case the phase — arg(Zi(r)) flows to a value arctan(a//3) as determined by the expressions 
in Eq. (I140p . The important difference here is the fact that in this case we have a different 
charge configuration which should now respect the non-BPS constraints of Eq. (I143p . Clearly, 
similar results hold for the second and third qubits playing a special role. Notice also that 
the quantities \Za\ where a = 1,2,3 occurring in the expressions of the three-qubit states 
like Eqs. (I14ip and (I15ip are just the attractor values of the fake superpotential 

Wair) = \Za{T)\. (154) 

For the |-BPS case a similar role is played by the quantity W(r) = 1^(1")!. Since 

Z(r) = -72^111, Zi = 72^110, Z2(r) = v^^ioi, ^3 = ^2^011 (155) 

with the remaining amplitudes arising by complex conjugation (see Eq. (J30l) ) we see that the 
fake superpotential in the relevant cases is related to the magnitudes of those amphtudes 
which are not dying out as the corresponding BPS or non-BPS flow approaches the horizon. 



27 



C. The non-BPS Z ^0 case 



The general non-BPS Z ^ case^i is extremely different. The general attractor ffow 
solution is 

exp(-4[/(r)) = ho{T)h{T)h2iT)h{T) - h\ 

y \v 



where 



- (159) 

are charge-dependent constants with a and (3 given by Eq. (11401) . The real constants 
satisfying the constraint 

ai + 02 + "3 = 0, (160) 
account for the flat directionsiii^. The harmonic functions now defined as 

hj(T) = bj + i-h)-^r, (161) 

giving rise to the quantities 

Q = hbK + hoha + 2b, (162) 

= hK - hoha, (163) 

= hbh^ + hoha - 2b, (164) 

C4 = exp(-2f/) = ^/hohih2h3 - b^, (165) 

also making their presence in Eqs. (I156p . 

One can obtain the non-BPS seed solution^i investigated in Sec.|V]and|Vl]as a special case 
of the general non-BPS Z ^ solution with the parameters ?a = f?a = \J ~pbp^ ; = 1 = ^a, 
eta = 0. Here b = B/Gn is the undressed version of the B field of the seed soUution. 
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Now we can present the horizon-limit of with some calculational steps needed for its 
derivation are left for the Appendix 



1 



2 — sgn(z/) cosh(a3 + if) cosh(a2 + V^) cosh(ai + (p) 

r/o|000) + r^illlO) + r^allOl) + r/slOll) + + /€i|001) + /talOlO) + K3IIOO) , (166) 



where 

kq = -13, 
and let u = sgn(i/) e'^ 



rja = i (— /3 sinh(— + 2(y9) + a cosh(— + 2ip)) , (167) 
Ka = sgn(i/) (/3 cosh(aa + ip) — a sinh(aa + ip)) , (168) 



(yj = In 



(169) 



Let us now consider some special charge- configurations giving rise to non-BPS Z ^ 
solutions. For the DO — DA configuration only the charges go and p"" are switched on. We 
consider the case when qo < and p°- > 0. Then J4 = Aq^p^p^p^ < 0, a = 0, /3 = 2p^p'^p^ > 0, 
z/ = 1, and 



-Aq^p^p^p^^ 



1^) 

2A/cosh(a3) cosh(a2) cosh(a;i) 
?sinh(ai)|110) + z sinh(a2)|101) + i sinh(a3)|011) 

-|lll)+cosh(ai)|001) + cosh(a2)|010) + cosh(a3)|100) 
As a special case when cti = 02 = 03 = one gets 

|§) = (-4gopW)^-f|001) + |010) + |100) - |111) 



(170) 



(171) 



i.e. we get back to the state at the horizon obtained for the seed solution in Eq. ( I124p . 

Now we consider the dual case of the D2 — D6 charge-configurationr^'^i^. After choosing 
p° > 0, ga > 0, I4 = — 4p°gig2g3 < 0, a = ^/\h\p'^ > 0, /3 = 0, = —1, one obtains 



1^) 



(4p°gig2g3) 



2A/cosh(a3) cosh(a;2) cosh(ai) 
-i|000) +icosh(ai)|110) +icosh(a2)|101) +icosh(a3)|011) 

+ sinh(ai)|001) + sinh(a2)|010) + sinh(a3)|100) 



(172) 
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Specially, when ai = 02 = «3 = 

1^) = i(4pVg2g3)^^ [|iio) + |ioi) + |oii) - |ooo)] . (173) 

Comparing Eqs. fll7ip and (11731) one can see that for vanishing flat directions the DO — DA 
amplitudes are real and the D2 — DQ ones are purely imaginary. Moreover, these cases are 
dual in the sense that they are related by the bit flip operation cti cti ai . Notice also that 



the norms of these states give ^^/—Aqop^^^]^ and \/ Ap^qiq2Q'i apart from the fact whether the 
flat directions are vanishing or not. These quantities multiplied by tt give the macroscopic 
black hole entropy^. 

It is interesting to analyse the effect of the asymptotic data on such states on the horizon. 
More precisely we are interested in those changes that leave the entropy (i.e. the norm of the 
state) invariant. As an example let us consider the D2 — DQ case. In the hope to have an 
effect merely on the relative phases of the state at the horizon we can adjust the asymptotic 
values for the charges and the aa parameters (flat directions) coming from the moduli. Other 
information coming from the asymptotic moduli are swallowed by the attractor mechanism. 

Let us change the signs of the charges qi,q2 and ^3, in such a way that the constraint 
P^qiq2q3 > is not changed. Then one can show'' that the possibilities for |\I') are 

\^)m,m,m, = 2(4Aig2g3)^^ [milllO) + m2|101) + malOll) - |000)] . (174) 

where 

(m3,m2,mi) G {(+ + +),(+--),(- + -),(--+)}■ (175) 

Hence although these changes are not affecting the black hole entropy, they have an effect 
on the particular form of the state. As one can check the possible changes giving rise to the 
four states of Eq. (11741) can be represented by phase flip error operators as 

(T3 (g) aa (g) /, 0-3® J® (T3, /®cr3®(T3, (176) 

where for example 

|^)__+ = (a3®a3®/)|^)+++. (177) 

Alternatively we can apply the corresponding bit flip error operators containing cji on the 
Fourier transformed states. Notice that these four states are all invariant under as^a^^a^. 
The result of this is that the subspace spanned by these states is invariant under an arbitrary 
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number of phase flips (or bit flips for the Fourier transformed subspace.) We can thus 
conclude that the effect of the change of this particular type of asymptotic data on the state 
at the horizon is the appearance of phase or bit flip errors associated with an invariant 
subspace spanned by the states of Eq. fll74p . 

Let us now remain in the non-BPS Z ^ charge orbit and fix the signs of the D2 — D6 
charges, but this time let us change the asymptotic values for the parameters of the flat 
directions from = to 7^ 0, ai + 0:2 + as = 0. In this case we see that the uniform 
structure of Eq. f ll73p deteriorates via the schematic transformation rule (neglecting the 
normalization factor) 

|000) ^ |000), |110) ^ coshailllO) -isinhailOOl), e.t.c. (178) 

Let us denote this new state i.e. the one of Eq. fll72p by \'^) 030201- Then one can show that 

1^)030201 = (^3 ® ^2 ® ^i)|^)+++ (179) 

where 

1/1 \ 

Ea = , ai + a2 + as = 0. (180) 

A/coshaa W sinh aa cosh aa I 

Hence the changes on the state |\^)+++ originating from the flat directions have the obvious 
interpretation of errors of more general kind depending on continuously changing parame- 
ters. Notice that we would have obtained the same state after changing the sign of the term 
cosh aa in the lower right corner of the matrix in Eq. fll80p . Such matrices G GL(2, C) 
in the limit a^ — )■ result in the phase flip error operators (Xs acting on the corresponding 
qubit we have already discussed. In quantum information theory the GL(2, C) operators 
acting on the qubits are called transformations associated with stochastic local operations 
and classical communication (SLOCC)^. It is amusing to see that though the error opera- 
tors E^ act locally but the constraint ai + a2 + as = (in the type IIA duality frame coming 
from deformations preserving the overall volume of T^) refers to the fact that they are not 
independent. In quantum information theory such constraints usually refer to an agreement 
between the parties effected via the use of classical channels. Finally in closing this subsec- 
tion we note that the normalized part of the attractor state of Eq. ( I173p i.e. |\&)+++ is just 
the one which can be used to establish a very striking version of Bell's Theorem.— 
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D. The DO - D6 case 



As another special subcase of the non-BPS Z one, finally we consider the DO — D6 
solution. This charge- configuration can only appear in the non-BPS regime, because I4 = 
— {p^qoY < independent of the signs of the charges. Originally, the general non-BPS Z 7^ 
solution was produced by an SL(2,]R)®^ U-duality transformation of the DO — D6 one^i"^. 
Due to the ^a, 0a parametrisation of this transformation, (see Eqs. (11581) and (11591) ) we can 
not produce neither the identity transformation nor the transformations that bring us back 
to the DO — D6 case with different charges. Hence we can not simply write the DO — D6 
charges into the corresponding formulae for the horizon-limit of ^I'(t). However, for the 
calculation of ^^(r) on the horizon we can directly use the original DO — D6 solutions^i 
instead. 

exp(-4f/(r)) = ho{r)h{T)h2{r)h{r) - h\ (181) 
S'^W = (182) 
nr) = (183) 



with the notation introduced in Eqs. (I162l) - (ll65p . and 



3 



184) 



Note, that on the horizon x*^ — > and — ?■ z/^. With this moduli the calculation of \E' is 
much easier, than in the general case. Finally a straightforvard calculation yields the result 

|§) = ^(-14)3 [-z|000) + |111)1 . (185) 
V2 L 

Note that this state is independent of the Ua parameters of fiat directions. 



VIII. CONCLUSIONS 



In this paper we have shown that in the special case of the STU model the attractor 
mechanism for extremal, static and spherically symmetric BPS and non-BPS black hole 
solutions can be cast in a form of a distillation procedure of entangled three-qubit states of 
special kind on the horizon. Such states are belonging to the so called GHZ-class featuring 
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maximum tripartite entanglement^. In obtaining this result our main calculational tool was 
to organize the charges, the moduli fields and the warp factor in a suitably defined three- 
qubit state (Eqs. (^^ and dZnD). This state is just lying on the GL(2,C)^3 orbit of a 

three-qubit "charge state" (Eq. fl2^ ) fixing the charge orbit to which the particular solution 
belongs. To conform with the well-known fact that the duality group of the STU model 
is not GL(2, C)^^ but SL(2,M)^^ our state is also satisfying a reality condition (Eq. (l3T|l ). 
On this three-qubit state the fiat covariant derivatives are acting as phase (sign) or bit fiip 
errors (Eq. (El])) depending on whether we express the sate in the computational basis or 
in its discrete Fourier transformed version (Eq. (I55|) -( 15^ ). The black hole potential can be 
expressed as the norm of our state (Eq. (ITTI) ). 

For spherically symmetric solutions such states are also displaying an explicit dependence 
on the radial coordinate r = ^ referring to the distance from the horizon. By solving the 
equations of motion for the moduli fields and warp factor we end up with a fiow of BPS or 
non BPS type depending on the charge configuration. Using the explicit forms of these fiows 
that has recently become available in the literature we can study the distillation procedure 
in detail. In order to illustrate how this distillation becomes unfolded as we approach the 
horizon we have chosen the recently discovered non-BPS seed solution. 

For such solutions we observed that the charge, moduli, and warp factor dependent state 
of Eq. f l76|) satisfies a system of first order differential equations (Eqs. f ll02p and fllOQp ) 
featuring the fake superpotential (Eq. (llOSp ). This observation conforms with the recent 
results on the first order formalism relating supergravity fiows to geodesic motion on the 
moduli space of the 3D dimensionally reduced theory^ii^. In the light of this connection 
it would be nice to elaborate further on this point and establish an entanglement based 
understanding of some of these results. 

For the non-BPS seed solutions we managed to demonstrate how a standard GHZ state 
at the horizon emerges from a state characterizing the fiow at the asymptotically Minkowski 
region. The attractor mechanism in this picture simply amounts to the fact that three 
amplitudes out of the seven nonequal nonvanishing ones of our three-qubit state are dying out 
as we approach the horizon (see Eqs. f ll23p and f ll24p ). The remaining amplitudes have the 
same magnitudes related to the macroscopic black hole entropy. The relative phase factors 
of these amplitudes are merely signs refiecting the srtucture of the fake superpotentiaUii^. 

In this paper we also conducted a detailed study on the structure of "attractor states". 
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By this term we denote the particular states we obtain from our r dependent ones after 
performing the r — )■ oo hmit. We have shown that the |-BPS and non-EPS Z = cases 
are very similar. The attractor states are of canonical GHZ form with the relative phases 
related to the phases of the central charge Z or the phases of the quantities Za, a = 1,2,3 
which are just the flat covariant derivatives of Z (Eqs. fll4ip and (11511) ). We observed that in 
these cases and also in the case of the Z seed solutions the fake superpotential is related 
to those amplitudes of our r dependent states which are not dying out as we approach the 
horizon. 

The new subtlety arising in the non-BPS Z ^ case is the appearance of fiat directions. 
As it is known for the DO — DQ case the parameters labelling the fiat directions are related to 
deformations of the volumes of the three tori preserving the overall volume of (in the 
type IIA duality frame). We have found that for this charge configuration the fiat directions 
are not making there presence in the corresponding attractor state. (Though they do appear 
in the particular form of the attractor values of the moduli.) However, for the most general 
charge configuration fiat directions do appear in the attractor states. In the special cases 
of the DO — DA and D2 — D6 systems we have shown that the effect of the fiat directions 
is to deteriorate the uniform structure of the corresponding attractor states obtained by 
starting the flow not in any of the fiat directions. It is known^ that the effect of changing 
the signs of the D2 and DA charges asymptotically results in phase or bit flip errors on the 
attractor states. By virtue of this the presence of fiat directions adds an additional twist 
to this picture. In particular we have demonstrated that fiat directions can entertain the 
possible interpreatation as errors of more general type (i.e. ones depending on continuously 
changing parameters) acting on attractor states. 

Now we comment on the possible physical relevance of our three-qubit states. Obviously 
our compressing of the variuos ingredients of the STU model in a three-qubit state at this 
stage is merely a nice way of understanding the structure of BPS and non-BPS solutions in 
the STU model. Notice however, that the attractor states are always just the ones that are 
connected to the structure of the fake superpotential. Indeed, the fake superpotential in our 
examples turned out to be related to those amplitudes of our three-qubit states which are 
not dying out during the process of distillation. It is known that upon quantization of the 
radial evolution of the moduli^^ results in a semiclassical wave function the phase of which 
is featuring the quantity e^W which is just formed out of the aforementioned amplitudes of 
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our three-qubit state of Eq. (1761) . (See also Eqs. (11050 . (I107P and (llOSp .) This observation 
might give a clue towards getting a deeper insight into the physical meaning of our GHZ-like 
states. 

Finally notice that although being very special, the STU model captures the essential 
features also of extremal black holes in the = 4, 8 theories. Moreover many of its features 
generalize well to other black hole solutions (such as those arising from CY compactifica- 
tions). In this respect we just remark that the maximal N = 8, d = A supergravity has seven 
STU subsectors corresponding to its consistent truncations. In the corresponding extremal 
black hole solution context this observation has already been related to systems exhibiting 
tripartite entanglement of seven qubits^i^. It would be interesting to study distillation issues 
for this more general scenario using the ideas as developed in this paper. 



Appendix A: Calculating ^'(t) on the horizon (non-BPS Z ^ case) 

In this appendix we outline the main steps leading to the explicit expression of \1/ on the 
horizon. First recall the explicit form of the non-BPS Z ^ attractor flow of Eqs. (I156p - 
(I165p . Using the (I158P and (I159P forms of and Qa, we can alternatively write this flow as 

= + (Ai) 

2{^pbpc _ pUg^j 2{py'' - p^qa) 
rir) = , C^y{r). (A2) 

Here the r dependent terms are 

^y^"^ = ^j^cirkkrcj- ^^^^ 

Since the moduh can be written as x^i^r) = A" + B^CI^t) and y"(r) = B^C^ir), the 
transformation Sa of Eq. 1^5^ can be expressed as 

Sa = -^CaBaAa (A5) 



,ya 



where 

Ca= I 

-C" 1 




= I ° I , Ba=\ " 1 , = I ^ " I . (A6) 
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After these preliminaries the transformation 6*3 ® 5*2 ® 5*1 can be carried out in three 
steps. First in order to use the (11401) definition of a and (3 also in this non-BPS case, the 
identity f ll37p have to be changed as 

4(pV - p%,){py - /g2)(py - /gs) = (A7) 
due to /4(r) < 0. By virtue of this we can perform the first two transformation 

(S3^3®S2^2®SiA)|r) = 

1 -h 



a(|000) + |011) + |101) + |110)) + + |100) + |010) + |001)) 

^^(/3 + a)|000)-l(/3-a)|lll) 



^2_,.L2- 2- —J 

where |0) and |1) are the ( |39l) Hadamard transformed states. On this latter form the 
transformation (C3 C2 ® Ci) acts readily. The r-dependency appears only in Ca- In the 
horizon- limit we obtain 



lim C; 



lim 

T — inn " 



1 



(A9) 



with these formulae we get the result of f ll66l) - fll67p . 
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